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1. Introduction 

Random Matrix Theory, introduced to theoretical physics by Wigner and Dyson E] 
more than 60 years ago, has recently seen important new developments prompted by 
its application to several different problems. Of particular interest are generalizations 
to ensembles involving two independent matrices, such as 2 hermitian random matrix 
theory (2HRM) and normal random matrix theory (NRMT). Rigorous results obtained 
in 2HRM over the last few years jHl IH El El El EI have led to important 
relations concerning bi-orthogonal polynomials, the Riemann-Hilbert problem, and the 
Kadomtsev-Petviashvilii (KP) hierarchy of integrable differential equations. Results of 
similar nature have been obtained for normal matrix ensembles, sometimes with simple 
geometric interpretations fOl El El El El El CSj , relevant to conformal maps in two 
dimensions [HlElEj- 

Critical points of hermitian random matrix ensembles have been studied intensively 
because of their important relations to 2D quantum gravity and string theory [201 
1^ 122 I2ni 1211 12S|; a similar analysis for generalized models of two matrices is under 
development. In NRMT, the evolution of the system towards the critical point has been 
given a clear hydrodynamic interpretation fOl El El El El El Ej in terms of the 
celebrated Hele-Shaw free-moving boundary problem |2Sll2Ill2Hll2niinnilSIllS2llSn!. In 
this paper, we continue this analysis and explore the scaling behavior of the ensemble 
at a critical point. 
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2. Normal Random Matrix Ensembles 

Definitions We briefly recall the basic concepts of NRMT ESI EEl IE] : the object 
of study is the ensemble of N x N normal matrices M {i.e. which commute with their 
Hermitian conjugates, [M, M'''] = 0), with statistical weight given by 

eb^mM,M^)d^^M), (1) 

and we choose for the present work the function W to be of the form 

W = -\z\^ + V{z) + V{z), (2) 

where V{z) is a holomorphic function in a domain which includes the support of 
eigenvalues. In (0), is an area parameter, and the measure of integration over 
normal matrices is induced by the flat metric on the space of all complex matrices. 
Upon integrating out angular degrees of freedom, the joint probability distribution of 
eigenvalues zi, . . . , is expressed as 

N 

i^^\^N{z)\'l[^^'^^'"'^^^'zj (3) 



Here d'^Zj = dxjdi/j for Zj = Xj + iyj, A]\r{z) = det(z* '^)i<i,j<N = ni>j(-^i ~ ^j) the 
Vandermonde determinant, and 



N 

is a normalization factor, the partition function of the matrix model. 



rN = ^ [ \A,,iz)\^ He'^^^^^-^Uh, (4) 



Droplets of eigenvalues It has been known for a long time that in a proper large N 
limit {h — s> 0, to = fixed), the eigenvalues of M densely occupy a domain D in the 
complex plane. The first rigorous result was obtained by Ginibre [32] in the form of 
Circular Law, later generalized by Girko to Elliptic Law [HZ], illustrated in Figure ^ 
In this case, the function V{z) introduced in (|21) is a quadratic polynomial, and the 
corresponding matrix model is referred to as Gaussian. Deformations of the droplet 
may be introduced by adding higher order terms to the function V{z), provided they 
are small compared to the unperturbed potential, in a domain including the droplet [33] . 
In the region of validity, a power expansion of V{z) = ^j^yi tkz'' is expressed through 
the exterior harmonic moments of the droplet of eigenvalues (I15|) [TH I15j . 



Wave functions and complex orthogonal polynomials In this section we specify the 
potential ^ such that it properly defines a scalar product for analytic functions in the 
sense of Bargmann |H8| IHH] . The wave functions and orthogonal polynomials are defined 
through 

ijjn{z) = e'^-^Xniz), and Xn{z) = e^^^'^Pniz), (5) 
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Figure 1. The distribution of eigenvalues for the Gaussian potential. The droplet is 
an ellipse with quadrupole moment 2|t2| and area T^hN . 



where the holomorphic polynomials Pn{z) are orthogonal in the complex plane with 
weight e~['^l^"^''^^^(^)l/'^. They obey a set of differential equations with respect to the 
argument z, and recurrence relations with respect to the degree n. Multiplication by z 
may be represented in the basis of Xn through a semi-infinite lower triangular matrix 
with one adjacent upper diagonal, Lnm = for m > n + 1: 

LnmXm{z) = ZXn{z) (6) 

(summation over repeated indices is implied), and the differentiation dz is represented 
by an upper triangular matrix with one adjacent lower diagonal. Integrating by parts 
the matrix elements of dz-, we have: 

(^L^^YimXm ^9zXni (7) 

where is the Hermitian conjugate operator. 

Operators L, L\ may also be represented in the basis of the shift operator W 
defined through Wxn = Xn+i, and become 

L = KW + f/(^) W-^ = W~^R + W'^U'-''^ , (8) 

k>0 k>0 

where R, U^^^ are diagonal semi-infinite matrices of elements Rnm = ^n^nm-, u!^''^ — 



nm 



Un^Snm, n,m > 0, and 6nm is the Kronecker symbol. Acting on x„, we obtain the 
commutation relation (the "the string equation", compatibility of Eqs. © and ((Tj)) 

[L\L]=h. (9) 

The string equation provides relations between the coefficients r„ and uii!^\ In particular, 
its diagonal part reads 

nh = rl-Yi2\^i'lf- (10) 

k>l p=l 

Example: Normal Gaussian Ensemble For the case of Gaussian potential V{z) = t2z'^, 
2|t2| < 1, orthogonal polynomials P„ are given by complex Hermite polynomials of 
scaled variable j4]j|. Recurrence relations (jH)) become two-term: 

and equations (01101) S^^^ the parameters rn,Un recursively through 

2t2 = ^, nh = rl-\u^^,\^ (12) 
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Figure 2. Growth in the non-critical Ginibre-Girko ensemble. The shape parameter 
of the ellipse \t2 \ < 1/2, so the normalized area = nh may be increased indefinitely 
without encountering any singularity. 



where the second equation is the area formula (jlOj) . Defining the two-component vector 
X = (Xn-i, XnY, we may express the action of the shift operator on x as 



1 



The operator V' also acquires a reduced 2x2 matrix representation: 



Z— — Un+l — 



I ^ _ "n+2«n+l Ur,+2 " 1 Xn+l' ^^^^ 



— 

" r„+i r„+i 

With the help of (HID, >C„ reads 



z{2t2y' -r^{2t2)-\l - Mhl^) 
r„(l-4|t2n 2t2Z 



Complex curve Since the operator L''^ represents multiplication by z, it is possible to 
define for each value of ra a complex curve fn{z, 5) = by solving the eigenvalue equation 
det[Cn{z) — z] = and setting z = z; we obtain an ellipse of equation 

.^-fe.= + *V)^^-nfti^ = 0. (14) 

with the quadrupole moment 2\t2\ and area vrto = vrn^. This is precisely the boundary 
of the domain D„ filled by eigenvalues of the matrix model, at fixed h, t2 and n (or 
normalized area to)- 

The geometrical meaning of the complex curve ()14|) is straightforward: at fixed 
shape parameter t2 and area parameter increasing n yields elliptical domains that 
represent the support of the corresponding n^n Gaussian model. A remarkable feature 
of this process (represented in Figure |21 and labeled growth in our previous works 
|14 | ITH t ITH]) is that it preserves the external harmonic moments of the domain D^, 

tfc(n) =tfc(n-l), tk{n) = -\ [ k > 1. (15) 



Jc\o„ z^ 

The only harmonic moment which changes in this process is the normalized area 
to = ^ J (Pz, and it increases in increments of h (hence the meaning of h as quantum 
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Figure 3. Deformation of the Gaussian ensemble and formation of boundary 
singularities. 



of area). We may say that the growth of the NRM ensemble consists of increasing the 
area of the domain by multiples of while preserving all the other external harmonic 
moments. The continuum version of this process, known as Laplacian Growth, is a 
famous problem of complex analysis. It arises in the two-dimensional hydrodynamics of 
two non-mixing fluids, one inviscid and the other viscous, upon neglecting the effects of 
surface tension (the Hele-Shaw problem |2nil2Zll2Hll2niEniinilin21ISS!). 

Deformations and critical points of the ensemble One of the most important properties 
of Laplacian Growth is that, with the exception of special choices for the external 
harmonic moments tk, the growth ends at a finite (critical) value of the area to = tc, 
when cusp-like singularities form on the boundary of the domain. 

Laplacian Growth of a simply-connected domain D(to) can be described jTUllIIllIl] 
through a conformal map z{w, to) which takes C\]Bi into C\D(to), where Bi is the unit 
disk centered at the origin in the w— plane. Therefore, the conformal map has the form 

z{w, to) = r{to)w + ^ Uk{to)w~^, (16) 

fc>0 

and the coefficients are functions of time (normalized area) to-, such that 

r / , \ -/ -1 , \T <ie/ \ dz dz dz dz 
M^,to),z>-,to)} = ..^----- 

The Poisson bracket (fT7j) encodes the infinite number of conservation laws 



1. (17) 



dtk _ , , If d^z 

— - U, tk{to) - --T / 



as well as the classical area formula 



r'-Y,kW?- 



As mentioned in the previous paragraph, we may regard Laplacian Growth as the 
continuum limit of a corresponding NRM ensemble, sharing the same set of exterior 
harmonic moments. For instance, choosing the function V{z) of the form 

n^)=«log(^)=E^-^ ^e(«)<0, (19) 
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corresponds to harmonic moments tk = and the exterior of the droplet is given by 
the map ^1] 



For suitable values of a, the orthogonal polynomials corresponding to ()19|) are well 
defined, though they may obey complicated recurrence relations. In the continuum 
limit of the model, the droplet grows until its area reaches the critical value and a 
cusp forms on the boundary. Figure El Formation of critical points is best described 
using the complex curve associated with the conformal map ()20|). As indicated in |14j . 
it is a degenerate elliptic curve, with two branch points at zi^2 inside the domain and 
a double point outside. The critical point on the boundary appears when one of 
the branch points inside merges with the double point, leading to a cusp expressed in 
local coordinates as ~ ip'. Restoring finite values for n and h is equivalent to a 
discretization of the Laplacian Growth and lifts the degeneracy of the complex curve of 
the continuum limit ITT)] . 

Universality in the scaling region at critical points: a conjecture Detailed analysis 
of critical Hermitian ensembles indicates that the behavior of orthogonal polynomials 
in a specific region including the critical point (the scaling region)^ upon appropriate 
scaling of the degree n, is essentially independent of the bulk features of the ensemble. 
This universality property (a common working hypothesis in the physics of critical 
phenomena) is expected to occur for critical NRM ensembles as well - and is indeed 
easy to verify in critical Gaussian models, 2\t2\ = 1. Analytically, it means that by 
suitable scaling of the variables z,n: 



where is the location of the critical point and tc the critical area, the wave function 
will reveal a universal part Q which depends exclusively on the local singular 
geometry x^ ~ y'^ {p, q mutual primes) of the complex curve at the critical point. This 
conjecture is a subject of active research. Its main consequence is that in order to 
describe the scaling behavior for a certain choice of p, q, it is possible to replace a given 
ensemble with another which leads to the same type of critical point, though they may 
be very different at other length scales. More precisely, the scaling behavior of operators 
in the vicinity of singular points illustrated in Figure El is assumed to be identical to that 
of singular points in Figure IH (where V{z) = t^z^), although the two critical droplets are 
obtained starting from different potentials. A constructive argument for this method is 
under development ji^ . 




(20) 



where 




(21) 



n — s> oo, h ^ 0, nh = to, to = tc — h^jy, z = Zc + K^C, 
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Figure 4. Critical points obtained from the cubic potential V{z) — t^z^ . Area of the 
droplet is increased at fixed until it reaches a maximal, critical value. 




Figure 5. Critical boundary corresponding to V{z) — t^z^ . 



3. Scaling at critical points of normal matrix ensembles 

In the remainder of the paper we analyze the regularization of Laplacian Growth for a 
critical point of type p = 3, g = 2, by discretization of the conformal map as described 
in the previous section. For simplicity, we start from the conformal map corresponding 
to the potential V{z) = t^z^, which is the simplest model leading to the specified type of 
cusp. It should be noted that the analysis will be identical for any monomial potential 
V{z) = tnz"',n > 3; for every such map, n singular points of type p = 3, g = 2 will 
form simultaneously on the boundary. The critical boundary corresponding to n = 6 is 
shown in Figure El 

Painleve I as string equation We start from the Lax pair corresponding to the potential 
V{z)=t3Z^ 0: 

Li^n = rntpn+1 + Untpn-2, ^Vn = r^-li^n-l + Mn+2^n+2- (22) 

The string equation Q [L'^ , L] = h translates into 

+ (r„_3U„ - rn-lUn-l)iJn-3 = ^^n- (23) 

Identifying the coefficients gives 

{rl - \u\l+2 - - (r^i - \u\l+i - \u\l) = h, (24) 
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-313. 



(25) 



Equation (j211) gives the quantum area formula 

rl - {\u\l+2 + =nh, 
which together with the conservation law ()25|) leads to the discrete Painleve equation 



(26) 



[l-9\hnrl_, + rl 



+1 



nh. 



In the continuum limit, the equation becomes 



- 18|t3| V = to- 
The critical (maximal) area is given by 

1^ = 0, 36|t3|V,2 = 1. 

Choosing rc = 1 gives 6|t3| = 1 and tc = \. It also follows that 
r„_2r„_i ^ 3 
2 ' 2' 



Ur 



Introduce the notations 



Nh = tc, nh = to=tc + h^^'v, rt = l- K^^'uiv), z 



where a = \. We get (9„ = ffd^ and 



' n+k 



1 - h'^u - kh-'^uM h^^u 



(27) 
(28) 

(29) 

(30) 
(31) 

(32) 



where dot signifies derivative with respect to u. The scaling limit of the quantum area 
formula becomes 



giving at order h'^"' the Painleve equation 

u-2u^ = 4z/. 
Rescaling u C2U, v ^ civ gives the standard form 

u — 3m^ = z/, 
for C2 = 4c?, Scf = 3. 



(33) 

(34) 
(35) 
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Painleve as compatibility equation Use the modified wave functions (Pol represents the 
polynomial part) 

n-l 



and rewrite the equations for the Lax pair as 



L4)n = (pn+1 + 



n-2' n-l 



■>n-2 



^_l^r^-l 



+ 



Notice that using the shift operator W, the system can also be written 

Introduce the scaling -0 function through 

The action of Lax operators on ip gives the representation 
L 



Therefore, the action of C, is given by the sum of equations at order fv^"-. 
and the action of 9^ by their difference: 



^2 



Equivalently, wc can write 



h'X^-[{W + lf+{rl_,W-' + iy 



4, 



'C 



Expanding the shift operator in h leads to 
W 

and 



1 + Wd, + fl'"^ + ^2"^, 

^ D 



a; 



6 

Substituting into the equations for 9^ gives the system of equations 

6 b 6 



(36) 



(37) 
(38) 
(39) 
(40) 

(41) 

(42) 

(43) 
(44) 

(45) 
(46) 

(47) 



where primed variables are differentiated with respect to C,. The equations can be written 
in matrix form as 



(48) 
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A = 



u , 2(C+u){2C-u) 
6 ~'~ 9 



2C,-u 
6 





2(C+^) 



The compatibihty equations 

k-Q' = [Q, A] 
yield the Painleve equation derived in the previous section: 



and 



Thus, 



A 



[Q, A] 



'it' I 2(^u—iuu 
6" ~' 9 



Q' 




^ 

3 ^ 



2{C,+u)u _u 
9 6 



= A - g' - [g, A] 







6 9 3 








The only non-trivial element of the matrix gives 

u —Auii — 4 = 0, 
i.e. the Painleve equation derived in the previous section. 



1 
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(49) 



(50) 



(51) 



(52) 



(53) 



(54) 



Painleve and the non- degenerate spectral curve In the scaling limit, the spectral curve 
is defined by the eigenvalue equation for the operator 9^, 

d^^ = A(C)vI>, (55) 

or 

det[A(C) - A] = 0. (56) 
We can write it also explicitly as an elliptic curve. 



^2^/2C_m\Y2C 2m\ 2C u u'-2uu 
3 3/V3 3/ 3'6 36 



The critical points of the curve solve 
A+ = A_ = 0, 

or 



, 2z/ 



27 



0, 



(57) 



(5^ 



(59) 



where /i = 2^/3. Setting all derivatives to zero in and (fS^. we get the degenerate 
solutions 



Cl = -U = -V-2U, C2,3 



(60) 



Generic critical points of normal matrix ensembles 



11 




Figure 6. Physical solution for Painleve I, with no singularities along the negative 
real axis. 




Figure 7. Evolution of the non-degenerate curve (the leftmost boundary in each 
diagram) and degenerate curve (middle boundary) relative to the singular curve (right 
boundary), at (a) v — oo, (b) Vf) — v = 0(1) and (c) — v vo. 

We choose (up to exponential corrections) the solution to Painleve I which is free of 
poles along the negative real axis, Figure El and follow the evolution of the boundary as 
^0 ~^ tc, Figure H As one can see, the presence of derivative terms lifts the degeneracy 
of the curve, so that the boundary remains smooth even as the area reaches its critical 
value. 



(1) Quantum curve at u —>■ — oo In the case u — > —oo, the asymptotic expansion 
of the solution to Painleve equation reads 

u = V^--^ + 0{u~'), (61) 

subject to exponential corrections. The curve becomes (Figure [7| first diagram) 
nondegenerate, with simple critical points 



1 ^ 1 



- ^2,3 - + ± (62) 

(2) Quantum curve at u = Let u^uq) = define the smallest real solution of 
m(z/) = 0. Then the local expansion reads 

2 

u{iy) = (x{y - z/o) + 2v^{y - v^f + -{y - v^f + —{y - z/q)' + . . . (63) 

6 6 

Merging this regular expansion with the asymptote at v ^ — oo yields 

z/o = -2^ a = A + ^4= = -3 ■ 2i (64) 
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The discriminant of aX u = uq becomes = ^" "'g?''"'' > 0, so the equation has 
again one real solution and two complex conjugate roots. Moreover, since the free term 
is positive, it follows that the real solution is negative. The physical interpretation 
shows that the curve is smoothed out (Figure diagrams b and c) at m = 0, instead of 
forming the classical (2, 3) cusp given by the degenerate curve. 
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